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Abstract 

If extra spacetime dimensions and low-scale gravity exist, black holes will be produced in ob¬ 
servable collisions of elementary particles. For the next several years, ultra-high energy cosmic 
rays provide the most promising window on this phenomenon. In particular, cosmic neutrinos can 
produce black holes deep in the Earth’s atmosphere, leading to quasi-horizontal giant air showers. 
We determine the sensitivity of cosmic ray detectors to black hole production and compare the 
results to other probes of extra dimensions. With n > 4 extra dimensions, current bounds on 
deeply penetrating showers from AGASA already provide the most stringent bound on low-scale 
gravity, requiring a fundamental Planck scale Mjj > 1.3 — 1.8 TeV. The Auger Observatory will 
probe Md as large as 4 TeV and may observe on the order of a hundred black holes in 5 years. We 
also consider the implications of angular momentum and possible exponentially suppressed parton 
cross sections; including these effects, large black hole rates are still possible. Finally, we demon¬ 
strate that even if only a few black hole events are observed, a standard model interpretation may 
be excluded by comparison with Earth-skimming neutrino rates. 
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I. INTRODUCTION 


Tiny black holes (BHs) can be produced in particle collisions with center-of-mass energies 
above the fundamental scale of gravity |Lj, ||, where they should be well-described semi- 
classically and thermodynamically In conventional 4-dimensional theories, viz., where 
the Planck scale ~ 10 19 GeV is fundamental and the weak scale ~ 1 TeV is derived from it via 
some dynamical mechanism, the study of such BHs is far beyond the realm of experimental 
particle physics. Over the last few years, however, physicists have begun exploring an 
alternative approach to the longstanding gauge hierarchy problem, wherein the weak scale 
becomes the fundamental scale of nature and the Planck scale is derived from this, with the 
hierarchy in scales a consequence of large or warped extra dimensions [|. |SJ. If this is the 
case, the fundamental scale of gravity can be 0( TeV), and BH production and evaporation 
may be observed in collisions of elementary particles B0.BBE3- 

If gravity indeed becomes strong at the TeV scale, ultra-high energy cosmic rays provide a 
powerful opportunity to probe BH production at super-Planckian energies [|ll]] . Cosmic rays 
with energies ~ 10 19 eV have been observed |T2|. They interact in the Earth’s atmosphere 
and crust with center-of-mass energies ~ 100 TeV, far beyond the reach of present and 
planned man-made colliders. These cosmic rays may therefore produce BHs, allowing cosmic 
ray detectors to test the existence of TeV-scale gravity and extra dimensions by searching for 
evidence of BH production |ll], [13]. [h|, [L5|, 0. A particularly promising signal is provided 


by ultra-high energy cosmic neutrinos, which may produce BHs with cross sections two or 
more orders of magnitude above their standard model (SM) interactions. These BHs will 
decay promptly in a thermal distribution of SM particles. Of the order of a hundred BH 
events may be detected at the Auger Observatory [UJ as quasi-horizontal, deeply penetrating 
showers with distinctive properties Jl3|]. The possibility of BH production by cosmic rays 


supplements possible sub-Planckian signatures of low-scale gravity JT7|, [18], [L£ 

In this article we extend previous work to derive bounds from the non-observation of 
BH-initiated showers in current data at the Akeno Giant Air Shower Array (AGASA). We 
also extend previous analyses of BH discovery prospects at Auger, and discuss in detail the 
possibility of distinguishing BH events from SM events. A preliminary version of some of 
these results was presented in Ref. [[Z^| . 

We begin in Sec. II with an overview of TeV-scale gravity. We collect and review existing 
bounds on the fundamental Planck scale in a uniform convention. In Sec. HI we discuss 
semiclassical BH production, including the effects of angular momentum and the production 
of Kerr BHs, as well as the proposed exponential suppression advocated by Voloshin |[23| |24f . 
This is followed in Secs. pA] and [V] by detailed discussions of cosmogenic neutrino fluxes and 
ground array experiments, respectively. 

Our results for event rates and new limits on the scale of higher-dimensional gravity are 
presented in Secs. 0 and Fill We begin with current data from AGASA. The AGASA 
Collaboration has already reported no significant signal for neutrino air showers during an 
observation time (live) of 1710.5 days ]^!J. Given the standard assumption of a geometric 


black hole cross section, we find that this data implies the most stringent bound on the 
fundamental Planck scale to date for n > 4 extra dimensions, exceeding limits derived |16 


from Fly’s Eye data ]26[ and also more stringent than the constraints from graviton emission 
and exchange obtained by the LEP [27j and D0 [^] Collaborations. In Sec. |VH| we then 
consider the prospects for BH production at the Auger Observatory. Tens of black holes 
may be observed per year; conversely, non-observation of BHs will imply bounds as large as 
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4 TeV on the fundamental Planck scale. 

In Sec. |Vffl| we note that comparison to Earth-skimming neutrino event rates 
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31, 32 


in Ref. ITl 


allows one to distinguish BH events from SM events. This point was noted already 
but was not considered in Ref. |16|], leading to weaker conclusions. Here, 


we consider this point quantitatively and find that, even with a handful of BH events, a 
SM explanation may be excluded based on event rates alone. If seen, black holes created 
by cosmic rays will provide the first evidence for extra dimensions and TeV-scale gravity, 
initiating an era of detailed study of black hole properties at both cosmic ray detectors and 

m ig, 


future colliders, such as the LHC 
Sec. m 


Our conclusions are collected in 


II. EXISTING LIMITS ON LOW-SCALE GRAVITY 


Depending on the dimensionality and the particular form of spacetime, the gauge hierar¬ 
chy problem may be re-expressed as a hierarchy in length scales. In the canonical example [[| , 
spacetime is a direct product of a non-compact 4-dimensional spacetime manifold and a flat 
spatial n-torus of common linear size 27T r c and volume V n = (27rr c ) n . Only gravity prop¬ 
agates in the full (4+n)-dimensional spacetime; all others fields are confined to a 3-brane 
extended in the non-compact dimensions. Here, the low energy 4-dimensional Planck scale 
Mpi is related to the fundamental scale of gravity in (4 +n) dimensions, M*, according to 


Ml, = M 2 + n V n = V n /G {4+n) , 


( 1 ) 


with G( 4+n ) defined by the (4 + n) dimensional Einstein field equation Rab ~ \9ab = 
— 87rG , (4 +n ) Tab- In what follows it will be convenient to work with the mass scale |38 


M d = [(2vr) ri /87r] 1/ ( n+2) M* , D = 4 + n . 


( 2 ) 


If r c is significantly larger than the Planck length, a hierarchy is introduced between Mp\ 
and Mb, and gravity becomes strong in the entire (4 + n)-dimensional spacetime at the 
scale Mb far below the conventional Planck scale Mpi ~ 10 19 GeV. Our conclusions will 
be essentially unchanged for more general “asymmetric” compactifications, with, e.g., p 
“small” dimensions with sizes < TeV -1 and n e g — n—p large extra dimensions |5D]. (Note, 
however, that in this case, the production of brane configurations wrapped around small 
extra dimensions may be competitive with black hole production J4(J.) Many of our results 
for black hole production and detection also apply for warped compactifications || in which 
the curvature length is much larger than a TeV' 1 . Hereafter, we will focus our discussion on 
bounds in flat compactihcation scenarios. In the figures, for n — 1 results, warped scenarios 
are implicit. 


A. Bounds from Newtonian gravity 

The provocative new features of these scenarios have motivated many phenomenological 
studies to assess their experimental viability. Naturally, the most obvious consequence of the 
existence of large compact dimensions is the deviation from Newtonian gravity at distances 
of order r c . For n — 1 and Mb ~ 1 TeV, r c ~ 10 13 cm, implying deviations from Newtonian 
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gravity over solar system distances, so this case is empirically excluded. For n — 2, sub¬ 
millimeter tests of the gravitational inverse-square law constrain Mo > 1.6 TeV [411. For 
n > 3, r c becomes microscopic and therefore eludes the search for deviations in gravitational 
measurements. 


B. Astrophysical bounds 


In the presence of large compact dimensions, however, the effects of gravity are enhanced 
at high energies, due to the accessibility of numerous excited states of the graviton (referred 
to as Kaluza-Klein (KK) gravitons j42|), corresponding to excitations of the graviton held 
in the compactihed dimensions. For low numbers of extra dimensions, by far the most 
restrictive limits on the radii of large compact dimensions come from the effects of KK 
graviton emission on cooling of supernovae, and from neutron star heating by KK decays |43 . 
For n = 2 the latter requires Mo > 600 — 1800 TeV, far above the weak scale; for n — 3, 
the bound is Mo > 10 — 100 TeV. These limits apply only for the situation where all extra 
dimensions have the same compactification radius. In the general case, the bounds could be 
less restrictive. 


C. Collider bounds 


For n > 4 extra dimensions, only high energy collisions are useful as probes. The effects 
of direct graviton emission, including production of single photons or Z’s, were sought at 
101 . 


LEP 


The resulting bounds are fairly model-independent, as the relatively low energies 


at LEP imply a negligible dependence on the soft-brane damping factor discussed below. 
For n = 4 (6), these null results imply M D > 870 (610) GeV [J27J]. 

The effects of low-scale gravity can also be seen through virtual graviton effects. These 
are most stringently bounded by the D0 Collaboration, which recently reported [^] the Erst 
results for virtual graviton effects at a hadron collider. The data collected at y/s = 1.8 TeV 
for dielectron and diphoton production at the Tevatron agree well with the SM predictions 
and provide the most restrictive limits on an effective extra-dimensional Planck scale for 
n > 4. This scale (called Ay in |38|, and related to Ms in |45| , [161 ) simply parameterizes the 


KK graviton exchange amplitudes for these processes: except for the different conventions 
used, they simply convey the experimental limit in terms of an energy-independent four-point 
function. In the context of low-scale gravity, the effective scale depends on both G( 4+n ) and 
on an ultraviolet cutoff on the contributing KK modes |38] [45], | 


This cutoff represents 


the energy where emission of graviton modes from the brane into the extra dimensions are 
damped by the effects of a non-rigid brane, and it is expected to be of order 


In this work we will use a Gaussian cutoff [[47], |48|], which emerges if one includes in 
the interaction the brane Goldstone modes. With this cutoff, a form factor e - m2 / 2A2 i s 
introduced at each graviton-matter vertex, where m is the mass of the graviton and A 
parameterizes the cutoff. In real graviton emission processes, the effect of the cutoff is 
somewhat alleviated because of finite cuts on the missing energy. However, at LHC energies, 
the corrections become significant |F| in the expected region A < Mo, and are of order 100% 
when A/Mo — 0.5. For virtual processes, the s-channel diphoton or dielectron amplitude 
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( 3 ) 


has the form 38, 45 


A 

S(s) 

T 


S(s)T 

8n . ^ 



Tfjsr' 


_j_ r T~l [l r T~IV 

71 + 2 » u 


where the sum on £ denotes a sum over the KK graviton modes, labeled by an n-dimensional 
lattice vector £, with graviton masses m = \£\/r c . The first and second 7%/s are the stress 
tensors for the incoming qq, gg and outgoing e + e _ , 77 states, respectively, and is the par- 
ton center-of-mass energy. The sum on £ may be approximated by a continuous integration 
over KK masses, modified by the cutoff, with the result |38], [45] 


S(S) 


S n - 1 r°° m n 1 din e m2 / A2 

Mj +n Jo m 2 -s 


( 4 ) 


where explicit integration over the n — 1 angular variables leads to the factor S n _ 1 = 
27r n / 2 /T(n/2). The connection to an effective four-point contact interaction in Refs. |38| , 
|7o| . m\ i s made by setting .s — 0 in Eq. ([|). This allows an explicit evaluation of the 
integration over m, with the result (for n > 3) 


5(0) 


n /2 ^ 

7r ' - 

n — 2 



A \ n " 2 1 

~Md) 


( 5 ) 


In the last line we have used the convention of Ref. [38] to parameterize the four-point 


amplitude. At 95% CL, the Tevatron data require Ay > 1.2 TeV. With the use of Eq. 
this allows us to generate Table |. which shows the bounds on Mr, for n = 3,... ,7 and 
0.5 < A /Mb < 1 . It is important to note that (except for small variations for the case of 
Ref. [|46|], which permits a sign ambiguity in the amplitude) Table [J is independent of the 
conventions in /^, [%|, ^/. We can see that the lower bounds on M D depend on both n 
and A. Typically, M D m[n < 1 TeV. 


III. BH PRODUCTION IN PARTICLE COLLISIONS 

The preceding section discussed some potentially observable consequences of scenarios 
with TeV-scale gravity, and the limits on the scale of higher-dimensional gravity resulting 
from their non-observation. In particular, for n — 4 to 7, the quoted lower limit on Mr, 
comes from the non-observation at D0 of processes involving KK gravitons. The spectrum 
and interactions of KK gravitons are model-dependent, to an increasing degree at increasing 
scales above Mr,. Here we describe a more universal and model-independent prediction of 
low-scale gravity scenarios: the production in particle collisions of microscopic BHs. 
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TABLE I: 95% CL lower limits on Md from the D0 Collaboration at the Tevatron. 


A /M D 

Md ,min (TeV) 


n = 3 

n = 4 

n = 5 

n = 6 

71 = 7 

0.5 

0.98 

0.80 

0.70 

0.63 

0.58 

0.6 

1.02 

0.88 

0.80 

0.76 

0.73 

0.7 

1.06 

0.95 

0.90 

0.89 

0.88 

0.8 

1.10 

1.01 

1.00 

1.01 

1.04 

0.9 

1.13 

1.07 

1.09 

1.14 

1.21 

1.0 

1.16 

1.13 

1.18 

1.26 

1.38 


A. Geometric cross section 


It has been argued [§[] that BH formation should occur in the scattering of two incident 
particles when their impact parameter is approximately less than the Schwarzschild radius 
of a BH of mass equal to their center-of-mass energy V§- This suggests a geometric cross 
section 

( 6 ) 


a 


nr. 


where 




1 

" A/bh 

1 

1 + 71 

'2 n V n - 3)/ 2 r( w + 3 )' 

m d 

Im d \ 


71 + 2 


1 

1 + 71 


(7) 


is the radius of a Schwarzschild BH of mass Mbh = Vs 50, 51] in 4+n dimensions. Even if 


the incident particles are stuck on the SM brane, the black hole formed should be treated 
as a fully 4 +n dimensional object in an asymptotically Minkowskian spacetime, as long as 
r s is small compared to r c . 

The cross section Eq. (^) grows like s 1 /( n+1 ) ; more rapidly than any SM cross section. 
Thus, at energies sufficiently far above M^, BH production is expected to become the 
dominant process. In pp collisions at the LHC, rates as high as 10 8 events per year have 
been predicted in scenarios with Md ~ 1 TeV 10- 

In our investigation of BH production by cosmic rays, we will be most interested in 
collisions of neutrinos with atmospheric nucleons. Since, at the energy scale of interest, 
gravitational cross sections will be far smaller than the geometric area of a parton, we write 
the uN cross section as fO 


a(uN —> BH) = ^ 




dx fri(Vxs) fi(x, Q) , 


( 8 ) 


where s = 2rriNE ul the sum is over all partons in the nucleon, and the /) are parton distri¬ 
bution functions. We set the momentum transfer Q = min {Mbh, 10 TeV}, where the upper 
limit is from the CTEQ5M1 distribution functions |52|]. The cross section a(vN —> BH) 
is highly insensitive to the details of this choice 0- For example, choosing instead 
Q = minjrj 1 ,10 TeV} changes BH production rates by only 10% to 20%. For the conser¬ 
vative v fluxes considered below, our results are also insensitive to low x. (For concreteness, 
however, we extrapolate to x < 10” 5 assuming a power law behavior fi(x, Q ) oc ^'h+MQ)].) 
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FIG. 1: Cross sections cr{vN —> BH) for n = 1,... , 7 from below for Md = 1 TeV, x m i n = 1 (solid) 
and 3 (dashed), and parton cross sections (left) and -k (right). The SM cross section 
a(vN —> IX) (dotted) is also shown. 


Finally, is the minimal BH mass for which Eq. (|6]) is expected to be valid. The ap¬ 

propriate choice of Mgjf is subject to theoretical uncertainties, as discussed below. We 
define 


•£min 


M^/M d 


(9) 


and present results for various 1 < x min < 5. The dependence of our event rates on x m in is 
found to be rather mild. 

Cross sections for BH production by cosmic neutrinos are given in Fig. |l] for Mu — 1 TeV; 
they scale as 

2 +n 
1 + 71 

( 10 ) 


a(vN BH) cc 


Mb 


The SM cross section for uN —> IX is also included for comparison. (Note that cross sections 
rise with increasing n for fixed Mp, whereas they decrease for increasing n for fixed M*.) 
As noted in Ref. 0 , in contrast to the SM process, BH production is not suppressed by 


perturbative couplings and is enhanced by the sum over all partons, particularly the gluon. 
As a result of these effects, BH production may exceed deep inelastic scattering rates in the 
SM by two or more orders of magnitude. 

Although greatly reduced by the cross section for BH production, neutrino interaction 
lengths 

- '?) 


L — 1.7 x 10 krnwe 


are still far larger than the Earth’s atmospheric depth, which is only 0.36 kmwe even when 
traversed horizontally. Neutrinos therefore produce BHs uniformly at all atmospheric depths. 
As a result, the most promising signal of BH creation by cosmic rays will be quasi-horizontal 
showers initiated by neutrinos deep in the atmosphere. For showers with large enough zenith 
angles, the likelihood of interaction is maximized and the background from hadronic cosmic 
rays is eliminated, since these shower high in the atmosphere. 

Once produced, the BH will Hawking evaporate, provided the semiclassical approximation 
is valid. In this case, a Schwarzschild BH will behave like a thermodynamic system with 
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temperature 


and entropy 


T h = 


71+1 
An r. 


S = 


2 7 r ( n + 3 )/ 2 
' n+ 3 ' 


y.n+2 __ 


A n M bh r s 
n + 2 


( 12 ) 


(13) 


4G D r(2±i) 

According to the semiclassical description, a BH produced in a scattering event should 
be regarded as an intermediate state, which decays on a time scale 


1 


r 


Mi 


BH 


Md V Mo 


3+re 

l+n 


(14) 


and statistical fluctuations in the 
For M bh /M d = 5, Eq. ( [hip gives S 


Since r < 10 ~ 25 s for M D > 1 TeV and M B h ~ 10 TeV, the decay is effectively instantaneous. 
In the decay process, particles will be radiated into all available SM channels, into quanta 
with energies typically of order T B or above. These decays are predicted to lead to highly 
distinctive signals in collider events with high multiplicity, large transverse energy, 

hard leptons and jets, and a characteristic ratio of hadronic to leptonic activity. 

The magnitude of the entropy determines the validity of this picture. Thermal fluctua¬ 
tions due to particle emission are small when S' ^ 1 H 
microcanonical ensemble are small for y/S > 1 
ranging from 29 for n — A to 25 for n = 7. For /M D = 3 (or 1), S is about 13 (or 3) 
for n between 4 and 7. 

In searches for BH mediated events at colliders, it is essential to set X min high enough 
that the decay branching ratios predicted by the semiclassical picture of BH evaporation 
are reliable, as there are very large QCD backgrounds, and the extraction of signal from 
background relies on knowing the BH decay branching ratios reliably. This is especially true 
if one is attempting to determine discovery limits, where the overall rates for BH production 
are not necessarily large. Thus, in collider searches, a cutoff of x min = 5 or more may be 
appropriate. 

By contrast, the search for deeply penetrating quasi-horizontal showers initiated by BH 
decays can afford to be much less concerned with the details of the final state, since the 
background is, relative to colliders, almost nonexistent. As a result, the signal relies only 
on the existence of visible decay products, which, in this context, includes all particles 
other than neutrinos, muons, and gravitons. Indeed, there is very little about the final 
state, other than its total energy and to some degree its multiplicity and electromagnetic 
component HI , that we can reasonably expect to observe, since detailed reconstruction of 
prompt decay particles is not possible at cosmic ray detectors. Thus, it seems reasonable 
to choose a significantly lower value of M^ n than is needed for collider searches; in our 
estimates of rates below we will take x min as low as 1. While BHs of mass around M D 
will be outside the semiclassical regime, it seems quite reasonable to expect that they, 
or their stringy progenitors, will nevertheless decay visibly, whatever stringy or quantum 
gravitational description applies. 

As an illustration, we examine the scattering in the string regime. For n — 6 large 
extra dimensions, M s ~ g 4 J 4 M P i (M s = string scale, g s = string coupling). As shown 
in 0, the string cross section er saturates to ~ 1/M 2 for y/§ > M s /g s , (or in terms of Mpi, 
a ~ Mpj 2 g for y/§ > M Pi /g ^ 4 . As noted in []36| this matches onto the classical black hole 
cross section at an energy M s /g 2 s r\j Mpi/g 7 J 4 . Thus, if g s is not too small (implying a small 









hierarchy between M s and Mp\), the transition to the geometric cross section is rapid. In 
this work we adopt a minimum energy ~ 1 — 3 Mr ~ 3 — 9Mp\ (for n = 6), so that while 
there are probably stringy corrections, the cross section should be substantially geometric. 


B. Uncertainties in the Cross Section 

Although the details of the process by which BHs decay are not of great concern to us, 
the production process is of central importance, since our rates (and the lower limits we 
will be able to set on Mr) will depend directly on the form of the cross section. It should 
be emphasized that the heuristic arguments on which Eq. (||) is based only determine cr 
up to an overall factor of order one. Even at the classical level, our conclusions could be 
significantly affected by theoretical uncertainties in this factor, four sources of which we now 
discuss. 


1. Mass ejection 


Classical general relativity calculations indicate that the mass of a BH formed in a 
head-on collision is somewhat less (about 16% less) than the total center-of-mass energy. 
At least in four dimensions, this suggests that the formula Eq. (((]) should be modified 
by replacing r s (M BH ) by r s (0.84M BH ), leading to a slight reduction of a. Very recently, 
corresponding calculations in more than four dimensions have been presented [[55] . 


2. Angular momentum 


The analytic techniques used to study head-on collisions are not applicable to collisions 
at nonzero impact parameter. Thus the claim that a BH will be produced when b < r s (\/§) 
can only be expected to be true up to a numerical factor. 

One issue that arises at nonzero impact parameter that we can address is that the BHs 
formed will have angular momentum N34|. In particular, the Schwarzschild radius appearing 


in the formula Eq. ([/]) should more accurately be replaced by the radius of a Kerr BH of 
the appropriate angular momentum. This will alter the critical impact parameter at which 
a BH will form, for given s. For two particles each of energy E in the center-of-mass frame 
colliding with impact parameter b, the total angular momentum with respect to the center 
of mass is J = 2{b/2)E = 6M B h/2. So the maximum impact parameter at which a BH will 
form should occur at a value of b for which the radius 7y.(M, J = 5M B h/ 2) of a Kerr BH and 
b are equal. The Kerr radius satisfies [|50| 


M bb = c 


Ml— 1 


m k 


r{ + (n + 2) 2 J 2 /4Mj 


^BH 


= c n r. 


„n—l 


r 2 + (n + 2) 2 6 2 /16 


with c n an n-dependent constant. Setting ry = b we get 

1 + (n + 2) 2 /16 


M rh = c n r? +1 


'n' k 


(15) 


(16) 


Since for a Schwarzschild BH, M BB = c n r” +1 , the cross section in Eq. (|) should be corrected 
to 


a « 7rrfc(M BH , J) — 1 + (n + 2) 2 /16 


n+1 


7rr 2 (M BH ) . 


(17) 
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For 1 < n < 7, the correction factor is remarkably stable, ranging from 0.62 to 0.64. This 


result has been recently confirmed [55] in a classical analysis of black hole formation for 


collisions with non-zero impact parameter. We see, then, that including the effect of angular 
momentum also leads to a small reduction of er. 


3. Sub-relativistic limit 

While the corrections related to mass ejection and angular momentum both appear to 
decrease cr by factors of order 1, another potential correction to the naive cross section 
Eq. (p) could enhance it by a more-than-compensating factor. Namely, the critical impact 
parameter may be somewhat larger than the radius of the BH formed. An argument that 
supports this conjecture may be given in the case where the incident particles are sub- 
relativistic with rest mass approximately E = Mbh/2. In this case, the incident particles 
may be treated as BHs with mass M B h/ 2. If b < 2r< j (M BH /2), they will touch as they pass, 
and thus merge into a BH of mass M BB . In this regime, we expect as an approximate lower 
bound on a of 

cr > irb 2 = 7r[2r s (M BH /2)] 2 = 4 n/(n+1) 7rr s (M BH ) 2 • (18) 

For large n, the correction factor approaches 4. Of course, the situation may change con¬ 
siderably in the ultra-relativistic limit, but this estimate, in a limit that we understand, at 
least makes it plausible that the correct coefficient in Eq. (6i) might be larger than 1. Thus 
our choice to take a = tt r 2 s may well turn out to be conservative. 


4- Gravitational infall and capture 


It could also be argued that the cross section of Eq. (p|) is too small, because it supposes 
that a black hole only forms when the two particles come within a Schwarzschild radius of 
each other, when in fact we expect that gravitational collapse will occur for somewhat larger 
impact parameters as well. Another problem with Eq. (p|) is that while the cross section 
is measured with respect to the flat geometry of the asymptotic region, the Schwarzschild 
radius is a property of the highly curved region close to the singularity [^6 . 


A better measure of the cross-sectional area associated with a black hole of given mass, 
which overcomes these objections, is given by the classical cross section for photon cap¬ 
ture 0. if a beam of parallel light rays is sent in towards a Schwarzschild black hole from 


the asymptotically Minkowskian region of spacetime, the black hole’s classical cross section 
is defined to be the cross-sectional area of the portion of the beam that gets captured. In 
four dimensions, one finds from the geodesic equation (see, e.g., Ref. |57|) that 


a = 7 :b 2 c = 27t rG^Ml H 


(19) 


independent of the energies of the incoming photons. (The relevance of this cross section for 
black hole production has been independently argued in Ref. |58|].) The maximum impact 
parameter b c at which capture occurs is about 2.6 times as large as the Schwarzschild radius, 
and the cross section is enhanced by a factor of 27/4. A straightforward extension of this 
calculation to Schwarzschild black holes in 4 + n dimensions gives |7| 


b c 


. x 1. ( 3+n\ 

(3 + nRA+iJ 

-1 — r s ■ 

2 1 + n \Jl + n 


( 20 ) 
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For n = 1, the cross section is 4 times larger than the cross section of Eq. (^[), and for n — 7, 
it is still 87% larger. Thus, the enhancement suggested by this definition of the cross section 
could be significant and could easily offset other possible reductions. 1 


C. Exponential suppression? 


Quantum mechanical corrections to the amplitude for BH production may be even more 
significant than classical uncertainties. In particular, Voloshin f2d| . p!4| | has proposed that 
the cross section of Eq. (0) should be modified by an exponential suppression factor 


2 - I F 

a ~ 7T r s e b , 


( 21 ) 


with Ie the Euclidean Gibbons-Hawking action for the BH, which, in terms of the entropy 
of Eq. ( |T3| ) , is 

Ie = ——7 • (22) 


n + 1 

In part, Voloshin’s critique is based on previous attempts to calculate amplitudes for 
the production of classical field configurations (in which the multiplicity is greater than the 
inverse coupling) from initial quantum states. The intrinsically nonperturbative nature of 
such processes suggests employing an instanton-like approximation, which could lead to 
exponential suppression. Such an approach can be taken even for processes that are semi- 
classically allowed (such as multi-Higgs production pi M) that do not require tunneling 


in order to take place. The problem is not yet solved. For example, a recent lattice simu¬ 
lation f6l| shows no evidence for the enhancement of large multiplicity amplitudes manifest 
in perturbation theory, perhaps counter-indicating the formation of a classical field state in 
the quantum collision. 

We are cognizant of the uniqueness of gravitation (such as the onset of strong coupling 
for s > M^), and the support in favor of the geometric cross section based on classical 
calculations for both vanishing |54] and non-vanishing |55|] impact parameter. An additional 
supporting argument based on a string calculation has been given [|3(|] , and the applicability 


of CPT arguments when comparing black hole formation and decay, an element in Voloshin’s 
criticism |^3 


has been questioned 


Nevertheless, for completeness, we will also 


present results below for the exponentially suppressed cross section of Eq. For cosmic 
rays, we will see that even with Voloshin’s suppression factor included, useful bounds will 
emerge after 5 years of operation of Auger. 


IV. THE COSMOGENIC NEUTRINO FLUX 

Among the many possible sources of ultra-high energy neutrinos, the cosmogenic flux 
is the most reliable. This neutrino flux relies only on the assumption that the observed 


1 A somewhat more refined estimate of the cross section would take into account the rotation of the black 
hole. One may derive a geodesic equation for null geodesics in the equatorial plane of a rotating black hole 
(since the incoming particles are in the equatorial plane), and calculate the impact parameter 6 (Mbh, a) 
at infinite distance from a black hole of mass Mbh and rotation parameter a (as defined in Ref. |^(| |57j ). 
In four dimensions, the extremal value a = Mbh gives b = 2Mbh, reproducing the cross section of Eq. (|j). 
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extremely high energy cosmic rays contain nucleons and are primarily extragalactic in origin. 
If the charge of the primaries satisfies Z < 0( 1), as recently reported [0], an extragalactic 
origin is almost guaranteed, as the observed nearly isotropic angular distribution strongly 
disfavors galactic disk sources [Q. Moreover, even if the absence of the Greisen-Zatsepin- 
Kuzmin (GZK) cutoff |J3|] on cosmic ray energies is a reflection of our coincidental position 
near a nucleus/nucleon-emitting source, one still expects the full cosmogenic neutrino flux. 

Briefly, the argument for this is as follows m : the known astrophysical environments 
(within a few Mpc of the Earth) are not among the most powerful, but (in principle) can 
produce hadronic cosmic rays with the desired energies when parameters are stretched to 
their limits. Thus if these less powerful sources can accelerate particles above 10 20 eV, it must 
be that more powerful distant sources (like Fanaroff-Riley II radiogalaxies) can accelerate 
protons above photopion threshold, giving rise to the cosmogenic neutrino flux. Moreover, 
the approximately smooth power law behavior of the observed spectrum above 10 19 eV JT2 
seems to indicate that any “local source” contribution should be comparable to that of 
all other sources in the universe. Otherwise, one should invoke an apparently miraculous 
matching of spectra to account for the smoothness of the spectrum. This smoothness will 
provide the basis for obtaining the cosmogenic neutrino flux, as discussed in what follows. 

The chain reaction generating these cosmogenic neutrinos, triggered by GZK pion photo¬ 
production, is well known [[??], 0. The resulting neutrino flux depends critically on 
the cosmological evolution of the cosmic ray sources and on their proton injection spec¬ 
tra The high energy tail of the neutrino spectrum can also receive a significant 

contribution from semi-local sources, such as the Virgo cluster 0- Additionally, there is a 
weak dependence on the details of the cosmological expansion of the universe. For example, 
a small cosmological constant tends to increase the contribution to neutrino fluxes from 
higher redshifts |7T| . 

In our analysis we adopt the cosmogenic neutrino flux estimates of Protheroe and Johnson 
(PJ) |7D| . We consider their estimate with an injection spectrum with E cuto s = 3 x 

10 21 eV. In addition to and electron neutrinos also contribute to black hole production. 
In the high energy peak, the u e , and fluxes are nearly identical [0, and we include 
this v e flux in our analysis. The study of PJ incorporates the source cosmological evolution 


from estimates [0 


of the power per comoving volume injected in protons by powerful radio 

7| 


galaxies, taking into account the radio luminosity functions given in Ref, 
the resulting neutrino spectrum is shown in Fig. 0. The flux peaks around E 


The shape of 
i 2 x 10 17 eV, 


which is roughly the same energy suggested by other analyses following a source evolution 
proportional to (1 + z) 4 |69|, |7l||. To explore possible additional contributions from semi-local 


nucleon sources, we also consider below the cosmogenic neutrino flux estimates of Hill and 
Schramm (HS) 0], which are also given in Fig. A flux estimate of Stecker is also given 
there. As noted in 0], the PJ, HS, and Stecker fluxes all yield approximately the same 
rates for BH production. 

We stress that the PJ flux agrees with the most recent estimate 


7T] in the entire energy 


range, whereas the spectrum obtained in earlier calculations |)9| is somewhat narrower, 
probably as a result of different assumptions regarding the propagation of protons. The PJ 
analysis is performed within Friedmann cosmology with vanishing cosmological constant A, 
go = 0.5, and Ho = 75 km s^ 1 Mpc^ 1 , assuming an extragalactic magnetic field of 1 nG and 
a source spectrum proportional to TV 2 up to redshifts z = 9. The extension to cosmological 
models with A ^ 0 would not produce remarkable changes. For example, for Hm = 0.3 and 
12 a = 0.7, the neutrino flux is increased by a factor of < 1.7 
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FIG. 2: Cosmogenic v tl +u /t + u e fluxes from Protheroe and Johnson with energy cutoff of 3x 10 21 eV 
(solid) |7(J, Hill and Schramm (dashed) |72|], and previous estimate by Stecker without source 
evolution (dotted) |6^]. See text for discussion. 


V. ACCEPTANCE OF SURFACE ARRAYS FOR NEUTRINO SHOWERS 


Ultra-high energy cosmic neutrinos may be detected by ground arrays and fluorescence 
detectors on the surface of the Earth, as well as by space-based fluorescence detectors, and 
neutrino telescopes beneath the Earth’s surface. Here we concentrate on ground arrays, and 
consider two prominent examples: AGASA, the largest surface array currently in operation, 
and the Auger Observatory now under construction. 

AGASA consists of 111 scintillation detectors each of area 2.2 m 2 , spread over an area 
of 100 km 2 with 1 km spacing [|75|. The array detectors are connected and controlled 
through a sophisticated optical fiber network. The array also contains a number of shielded 
scintillation detectors which provide information about the muon content of the showers. 
The full AGASA experiment has been running since 1992, and has recorded the majority of 
events claimed to have energies above the GZK cutoff. 

The Auger Observatory is a hybrid experiment, with two sites (one in the northern 
hemisphere and one in the southern), each covering an area of 3000 km 2 and consisting of 
1600 particle detectors overviewed by 4 fluorescence detectors |76|. The surface array stations 
are cylindrical water Cerenkov detectors with area 10 m 2 , spaced 1.5 km from each other in an 
hexagonal grid. Event timing is made possible through global positioning system receivers. 
The optical system uses the fluorescence technique pioneered by the University of Utah’s 
Fly’s Eye detector [0- “Golden events,” events detected by both methods simultaneously, 
will be extremely valuable for experimental calibration. However, atmospheric fluorescence 
detection is possible only on clear, dark nights, and so the golden event rate is expected to 
be less than 10% of the total event rate. We consider only the ground array below. The full 
southern site is scheduled for completion in 2003. Its engineering array, at 1/40 of the full 
size, is now complete and is already detecting giant air showers |f78| . 

A surface array’s acceptance for neutrino detection may be expressed, in units of km 3 
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FIG. 3: Slant depths corresponding to various zenith angles 9. 


water equivalent steradians (km 3 we sr), as |79 


A(E) = S 


27r sin 9 dO 


Po 


Pwater 


■ e- h/H V{E,6,h) dh , 


(23) 


where S is the area of the ground array, po ~ 1.15 x 10~ 3 p wate r is the density of the at¬ 
mosphere at ground level, H cs 8 km, h max = 15 km, and V(E , 6 , h ) is the probability of 
detecting a shower with energy E and zenith angle 6 that begins at altitude h. The mini¬ 
mum zenith angle is set by the desire to separate deep neutrino-initiated showers from far 
showers initiated by hadronic primaries. Typically, a minimum zenith angle in the range 
60° < 6 m in < 75° is imposed. This range corresponds to atmospheric slant depths of 2000 to 
4000 g/cm 2 . (See Fig. |3|.) The maximum zenith angle 9 max varies from analysis to analysis. 
For example, in an analysis of fully contained showers, a value of 0 max below 90° is required 
for showers to deposit all of their energy within the array. 

Reliable Monte Carlo simulations to determine Auger’s acceptance for quasi-horizontal 
showers have been performed by several groups |79], |80|. Of course, the acceptance depends 
on the amount and type of energy generated by neutrino interactions in the atmosphere. 
For example, the charged current interaction —>■ fi + X produces a muon that carries 
approximately 80% of the incoming energy and is not detectable at Auger. Acceptances for 
both electromagnetic and hadronic showers have been determined in Ref. |3D|. BHs decay 
thermally, according to the number of degrees of freedom available, and so their decays are 
mainly hadronic. We therefore adopt the hadronic acceptance of Ref. [|80| including partially 
contained showers with zenith angles 6 > 75°. The acceptance is given in Fig. 0. Partially 
contained showers, where the shower axis does not pass through the array, do not contribute 
significantly to the Auger acceptance for shower energies below 10 10 GeV. 

The AGASA Collaboration has searched for deeply penetrating showers [^5|, |8l| . In these 
studies, they find that, for showers with energy above 10 10 GeV and the requisite zenith 


angle, the detection probability V becomes effectively 100% and independent of altitude [|25 
For these energies, then, Eq. (|23|) may be re-written as 


A(E > 10 10 GeV) « (Sfl) e ff(E) 


Po 


P water 


~ h ! H dh, 


(24) 
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FIG. 4: Ground array acceptances for quasi-horizontal air showers at the Auger Observatory (solid) 
and AG AS A (dashed). See text for discussion. 


where 


(sn) eS (E) = s 


2vr sin 6V{E,6) dd 


(25) 


is the “effective area x solid angle” 0- Acceptances for extremely energetic showers may 
then be quoted in terms of (ST2) e ff. The AGASA Collaboration has searched for deeply 
penetrating showers of any origin in Ref. [|8T|. They find none with energy above 10 10 GeV 


in 9.7 x 10 7 s of exposure. Given an upper bound of 2.44 events at 90% CL, then, they 


derive a flux limit for deeply penetrating showers of 1.9 x 10 10 km 2 sr 1 s l , which implies 
(sn) eS = 132 km 2 sr for quasi-horizontal air showers. Equivalently, given Eq. (|24|), the 
AGASA acceptance for neutrino initiated events is 


A(E > 10 10 GeV) « 1.0 km 3 we sr 


(26) 


This acceptance is roughly 30 times smaller than the neutrino acceptance of Auger, as 
one would naively guess from the ratio between the Auger and AGASA surface areas. For 
lower energies, since the separation between detectors is smaller at AGASA than at Auger, 
a conservative approach is to model the AGASA acceptance as that of Auger reduced by 
a factor of 30. We adopt this estimate for energies below 5 x 10 8 GeV, and interpolate 
smoothly between this and Eq. (^) for energies 5 x 10 s GeV < E u < 10 10 GeV. The 
resulting AGASA acceptance is shown in Fig. [4]. 


VI. NEW BOUNDS FROM AGASA 


Given the cross sections, apertures, and fluxes discussed above, the number of BHs de¬ 
tected by a given experiment is 


N = I dE v N a a{E v ) A{E V ) T , 


(27) 
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where A{E U ) is the experiment’s acceptance in cm 3 we sr, Na = 6.022 x 10 23 is Avogadro’s 
number, dQ/dE u is the source flux of neutrinos, and T is the running time of the detector. 
We now determine current bounds on BH production from the AGASA experiment. In the 
next section, we examine future prospects for BH detection at the Auger Observatory. 

The AGASA Collaboration has searched for deeply penetrating quasi-horizontal show¬ 
ers [^3]. The depth at which a shower is initiated is, of course, not directly measurable in 
a ground array. However, the electromagnetic components of far showers are extinguished 
by ground level, leaving only a muon component, whereas for deeply penetrating showers, 
both electromagnetic and muon components are detected. By exploiting this difference, 
deeply penetrating quasi-horizontal showers may be distinguished from showers induced by 
hadronic cosmic rays. 

Relative to showers with muon components only, showers with eletromagnetic components 
have charged particle densities that are more concentrated near the shower axis, and their 
shower fronts are more curved. The depth at shower maximum X max may then be determined 
through its correlation to two measurable quantities: rj, which parameterizes the lateral 
distribution of charged particles, and 6, which parameterizes the curvature of the shower 
front. The values of A may as determined by these correlations, denoted X ^ a x and A^ ax , are 
then required to be large to distinguish candidate neutrino events from showers induced by 
hadronic cosmic rays. 

In 1710.5 days of data recorded from December 1995 to November 2000, the AGASA 
Collaboration found 6 candidate events with X^ ax , A^ ax > 2500 g/cm 2 . The expected 
background from hadronic showers is , where the first uncertainty is from Monte 

Carlo statistics, and the second is systematic. Of the 6 candidate events, however, 5 have 
values of A^ iax and/or X^ iax that barely exceed 2500 g/cm 2 , and are well within AX max of 
this value, where AX max is the estimated precision with which X max can be reconstructed. 
The AGASA Collaboration thus concludes that there is no significant enhancement of deeply 
penetrating shower rates given the detector’s resolution. 

The AGASA results imply lower bounds on the scale of low-scale gravity, assuming the 
conservative cosmogenic fluxes of Sec. |IV|. For these fluxes, the expected rate for deeply 


penetrating showers at AGASA from SM neutrino interactions is about 0.02 events per 
year, and so negligible. Given 1 event that unambiguously passes all cuts, and the central 
value of 1.72 background events, the AGASA results imply an upper bound of 3.5 black hole 
events at 95% CL |83| . 

The 3.5 event contour is given for various dimensions n in Fig. [5|. For x m i n = 1, the 
absence of deeply penetrating showers in the AGASA data implies 


n = 4 : M d > 1.3 — 1.5 TeV 
n = 7 : M d > 1.6 — 1.8 TeV . 


(28) 


Results for x m i n = 3 are also given in Fig. ||. They imply Md > 1.0 — 1.1 TeV for n = 4, and 
Md > 1.1 —1.3 TeV for n = 7; even for s m i n = 3, these bounds are exceed or are competitive 
with all existing collider and astrophysical bounds. As argued in Sec. [HI A, 


x r 


= 1 


is a 


reasonable assumption for the present application, as the derivation of limits relies only on 
the assumption that BHs or their lighter progenitors with mass around Md decay visibly. 
This assumption is violated only if their decays are limited to neutrinos, gravitons, and 
muons. 

The range in Eq. (^) is from considering both PJ and HS fluxes. As noted in Ref. 


11 


the dependence of the bounds on variations in the evaluations of cosmogenic fluxes is weak. 
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FIG. 5: 95% CL lower bound on Mp from non-observation of quasi-horizontal air showers in 1710.5 
live days at AG AS A for x m i n = 1 (solid) and 3 (dashed), assuming the cosmogenic neutrino flux 
of Protheroe and Johnson (lower) and Hill and Schramm (upper). 


These bounds are conservative in that larger non-cosmogenic fluxes, as predicted by some 
models and as may be indicated by super-GZK cosmic rays, will strengthen them, possibly 
dramatically. Note also that we have neglected enhancements to cosmic neutrino interactions 
from sub-Planckian extra-dimensional physics, which are more model-dependent, but can 
only serve to strengthen these bounds. 

The bounds of Eq. ( |28|) are, of course, subject to the 0( 1) uncertainties inherent in 
the parton level cross section. Given this cross section, however, they are direct bounds 
on the fundamental Planck scale Mb, and are not subject to the uncertainties inherent in 
collider bounds, such as the choice of brane softening parameter A discussed in Sec. Eg Any 
comparison of collider and cosmic ray bounds is then subject to the independent uncertainties 
associated with each bound. Nevertheless, for n > 4, given the geometric BH cross section, 
the AGASA limit is more stringent that all existing collider bounds for all choices of A /Md < 
1 . 

Before leaving the AGASA results, we derive their implications for extra dimensions if 
taken at face value. Given 6 events with an expected background of 1.72 events, the expected 
signal is 0.86 to 11 events at 95% CL. The preferred region of the (n, Mjj) plane is given in 
Fig.g (In Fig. H and all following figures, we use the PJ flux. The HS flux yields slightly 
larger rates.) The evidence for BH production (or any other anomaly) is speculative, given 
the statistics and the peculiarities of the data noted above. However, this analysis shows 
the power of cosmic ray measurements for probing extra dimensions. The preferred Planck 
scales are not probed by any other experiment. At the same time, they will be thoroughly 
explored in the near future at larger cosmic ray experiments, such as the Auger Observatory, 
to which we now turn. 


17 



























FIG. 6: 95% CL upper and lower bounds on Md for various n, given 6 candidate events above 

a background of 1.72 in 1710.5 live days at AGASA, and ascribing the excess to BH production. 
We assume x m j n = 1 and the cosmogenic neutrino flux of Protheroe and Johnson. 

VII. FUTURE PROBES AT AUGER 


Given the apertures discussed in Sec. |V], it is a simple matter to estimate the BH event 
rate for Auger. The number of detected BH events are given in Fig. ^ for various n as a 
function of Md- The Auger ground array is expected to become fully operational in 2003. 
We assume a running time of 5 years, roughly the data expected before the LHC begins. 
For x min = 1, Auger will probe fundamental Planck scales as large as Md = 4 TeV. For 
M d ~ 1 TeV and n > 4, 100 BHs could be detected. 

Given the prospects for fairly high statistics, detailed BH studies are in principle possible. 
While BHs with mass near Mp are in some sense of the greatest interest, for detailed studies, 
one might first restrict attention to more massive BHs (more energetic showers), where the 
semi-classical description of BHs is expected to be justified. The distribution of BH masses 
in cosmic ray collisions is given in Fig. |8]. They are concentrated near Md, but the event 
rate is reduced by only 0( 1) factors for x min as large as 5. This contrasts strongly with the 
case at colliders, where there is little energy to spare, cross sections are suppressed by two 
parton distribution functions, and event rates are reduced by two orders of magnitude for 

Hi]. Total event rates for a: min = 3 are also given in Fig. 7|. 
100 BHs may be detected for Md near 1 TeV. 

The dependence of BH event rates on running time T is given in Fig. |Sj. The event rate 
contours rise rapidly at first — in even the first few months, Auger will be sensitive to values 
of Md beyond present experiments. 

If no enhancement of quasi-horizontal showers is seen, Auger will set stringent limits on 
low-scale gravity and scenarios with extra dimensions. To determine these limits, we again 


= 1 


x m in = 5 relative to x m ; n 
Even for x m i n = 3, we find that 


assume the cosmogenic fluxes of Sec. |TV| and that only SM sources of deeply penetrating 
showers are observed. In contrast to AGASA, SM neutrino interactions lead to observable 
rates - given the cross section of Fig. [lj. 0.5 events per year are expected. In addition, as at 
AGASA, hadronic showers may fake deeply penetrating showers. As noted above, the Auger 
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FIG. 7: Event rates in 5 years for the Auger ground array for .x m i n = 1 (left) and 3 (right). 
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FIG. 8: Event rates for BHs with mass above Mg^ n at the Auger ground array as a function of 
M^h'/Md for n = 1,... ,7 from below, assuming Mo = 1 TeV, 5 years running time, and parton 
cross section irr 2 s (solid) and irr 2 e~ lE (dashed). 


aperture of Sec. 0 assumes zenith angles 6 > 75°, corresponding to slant depths of X rnax > 
4000 g/cm 2 , significantly more stringent than for the AGASA study 0. Nevertheless, 
hadronic showers may be a significant background. We know of no detailed study, but 
consider the possibility of % background events from hadronic showers in 5 years below. 

Given these assumptions, the expected background in 5 years is roughly 2+ub events. To 
determine the expected limit on BH production we assume that 2 + ub deeply penetrating 
events are in fact observed. At 95% CL, then, the upper bound on signal events for ub = 0, 
5, and 10 is 4.7, 6.8, and 8.3 events, respectively, in Fig. |^. contours for these event rates are 
also given. We find that, for < 10, x min = 1, and n — 6, if no events above background 
are observed, Auger will extend current bounds on Mo to above 2 TeV after the first year 
of live time. After 5 years, for .x' min = 1, Auger will set a limit of Mo > 3 TeV for n > 4. In 
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FIG. 9: BH event rates at the Auger ground array for n = 6 and x' m ; n = 1 (left) and 3 (right). 
The dashed contours indicate the expected 95% CL lower bound on Mb in the absence of physics 
beyond the SM and assume ns = 0, 5, and 10 background events from hadronic showers (from 
above). The geometric cross section 7 xr 2 is assumed. 


conjunction with astrophysical bounds, this will require Mb > 3 TeV for all n, significantly 
straining attempts to identify the Planck scale with the weak scale in scenarios of large extra 
dimensions. Note that we have neglected model-dependent sub-Planckian effects that may 
increase the rates and strengthen the bounds presented here. 

Finally, we consider the impact of the proposed exponential suppression of BH production 
cross sections. In Fig. |^, we show the dependence on x m i n = M^ /Mb for black hole event 
rates including this suppression. For X min = 1, the exponential suppression is not particularly 
severe, reducing event rates by factors of 3 for large n. Of course, the impact is much larger 
for larger x min . In Fig. we show the number of BHs observed in time T for parton cross 
section 7ir 2 e~ lE . For x min = 1, Auger may still see tens of BHs in 5 years, and will extend 
current bounds to M D ~ 2.5 TeV. For x min = 3, the event rates are quite suppressed, but a 
few BH events are still observable in 5 years. 


VIII. DISTINGUISHING BLACK HOLES FROM SM EVENTS WITH EARTH- 
SKIMMING NEUTRINOS 


If an excess of quasi-horizontal showers is observed, how can it be identified as arising 
from BH events? After all, at first sight, an excess may arise simply from an enhanced flux. 
With sufficient statistics, a SM explanation may be excluded based on shower properties, as 
black hole showers differ markedly from those produced by SM charged and neutral current 
neutrino interactions |13|]. It may also be possible to confirm specific predictions of BH 
production by verifying Hawking radiation through correlations between A" max and shower 
energy 0 - 

It is also possible, however, to differentiate BH from SM events by considering additional 
constraints on ultra-high energy neutrino properties. In particular, comparison with Earth- 
skimming neutrino rates may allow one to distinguish BH and SM interpretations |LL|. In 
this section, we develop this possibility quantitatively, focusing on the question of excluding 
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FIG. 10: BH vent rates as in Fig. 
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but for parton cross section irr^e Ie 


a SM interpretation for BH events. We also comment briefly on the task of differentiating 
black hole events from other new physics possibilities at the end of this Section. 

At ultra-high energies, even the SM neutrino cross section is large enough that upward- 
going neutrinos are blocked by the Earth. However, neutrinos that skim the Earth, traveling 
at low angles along chords with lengths of order their interaction length, are not [29, |S(], 31 


32fl. These Earth-skimming neutrinos may then convert to charged leptons in the Earth’s 


crust, and the resulting charged leptons may emerge into the atmosphere, producing a 
signal in cosmic ray detectors. A schematic picture of such an event is given in the top 
panel of Fig. 0. The best signal is from r leptons. Unlike electrons that do not escape 
the Earth’s crust, or muons that do not produce any visible signal in the atmosphere, taus 
can travel for tens of km in rock, escape, and then decay in the atmosphere, leading to 
spectacular showers and observable rates of order 1 per year in both ground arrays 0 
and fluorescence detectors [|3Cj . The optimal angle for Earth-skimming neutrinos is energy- 
dependent. For E u ~ 10 8 (10 10 ) GeV, the optimal angle relative to the horizon is ~ 3° (1°). 
Given the angular resolution of cosmic ray detectors, these Earth-skimming events are easily 
differentiated from standard horizontal neutrino showers. 

The scenario changes radically in the presence of a significant cross section for BH pro¬ 
duction. First, BHs decay largely to hadrons, which do not escape the Earth. Such an event 
is pictured in the bottom panel of Fig. [□]. Of course, BHs also have a significant leptonic 
branching fraction, but leptons from BH decay carry only a fraction of the initial neutrino 
energy, and their detection rate is therefore highly suppressed. The probability of detecting 
BHs produced in the Earth by ground arrays and surface fluorescence detectors is therefore 
insignificant. Second, a sufficiently large BH cross section also depletes the original neutrino 
beam through absorption, leading to a substantial suppression of all Earth-skimming events, 
including those in the top panel of Fig. 0. 

To determine the effects of BH production on Earth-skimming rates, we consider here a 
simple analysis that is nevertheless sufficient to isolate the functional dependence of Earth- 
skimming rates on cross section parameters. The analysis extends the discussion of Ref. fl30| , 
where additional details and discussion may be found. 

Consider a flux of neutrinos with energy Eq. Given the high energies required for de- 
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FIG. 11: Top: A neutrino enters the Earth and converts into a charged lepton, which exits the 
Earth and may be detected. Bottom: A neutrino enters the Earth and produces a BH, which is 
captured in the Earth. 

tection, the most relevant energies are E 0 ~ 10 9 — 10 10 GeV, even for cosmogenic flux 
evaluations peaked at somewhat lower energies, and we may therefore limit the discussion 
to this rather narrow band of energy. Earth-skimming events occur in the Earth’s crust, and 
so the relevant neutrinos and taus sample only the Earth’s surface density, p s ps 2.65 g/crn 3 . 
In the SM, the neutrino’s path length is 

-^cc = [NaPs& cc] j (29) 

where Na — 6.022 x 10 23 g^ 1 and a^ c is the charged current cross sections for E v = E 0 . (We 
neglect neutral current interactions, which at these energies serve only to reduce the neutrino 
energy by approximately 20%.) For E 0 ~ 10 10 GeV, L£ c ~ 0(100) km. Supplemented by 
the possibility of BH production, the neutrino’s path length is 

^tot = \NaPs (^cc + ^bh)] j (30) 

where <Tb H is the BH production cross section for E u = E 0 . 

At these energies, the tail’s propagation length is determined not by its decay length but 
by its energy loss. The r lepton loses energy in the Earth according to 

dE 

—— = — (a T + (5 T E T )p s , (31) 

az 

where, for these energies, a T is negligible, and we take /3 r ~ 0.8 x 10 -6 cm 2 /g [|84|. The 
maximal path length for a detectable r is, then, 

L T = —- 111 (Emax/-E min ) , (32) 

PtPs 

where E max ps E 0 is the energy at which the tau is created, and E min is the minimal energy 
at which a r can be detected. For cosmogenic neutrino fluxes and other reasonable sources, 
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and the acceptances of typical cosmic ray detectors, taus cannot lose much energy and be 
detected. For E max /E min = 10, L T = 11 km. 

Given an isotropic v T + P T flux, the number of taus that emerge from the Earth with 
sufficient energy to be detected is proportional to an “effective solid angle” 


fl eff = J d cos 6 dcj) cos 9P(9,(j)) , 


where 


rt rf y 

P(0,<P) = / 6 [z- (l - IS)] 

Jo l;., 


(33) 


(34) 


is the probability for a neutrino with incident nadir angle 9 and azimuthal angle 0 to 
emerge as a detectable r. (In Eq. (|34|), for the reasons noted above, we have neglected the 
possibility of detectable signals from BH production by Earth-skimming neutrinos.) Here 
t = 2i?®cos0 is the chord length of the intersection of the neutrino’s trajectory with the 
Earth, with R® ps 6371 km the Earth’s radius. Evaluating the integrals, we find m 


fi eff = 2vr^ 

T V 

^CC 




2 Ras 



,-2 R®/E 


(35) 


At the relevant energies, the neutrino interaction length satisfies L” ot -C R®. In addition, 


for L^ ot L T , valid when the BH cross section is not very large, Eq. 


simplifies to 




eff 


ri/2 tt 

2vr - Ltot 


4 RlL 


(36) 


cc 


Equation (pq) gives the functional dependence of the Earth-skimming event rate on the 
BH cross section. This rate is, of course, also proportional to the source neutrino flux < F 1 ' 


at Eq. Finally, the constant of proportionality is determined by previous studies |29], |30 
where all the experimental issues entering tau detection have been included. Given these 
inputs, the number of Earth-skimming neutrino events detected in 5 years is 


Aes ~ Ces 


<jy 


cr: 


v 2 
CC 


cb u 

u'o 


(°cc + ^bh) 


(37) 


where Ces is the number of Earth-skimming events expected for the standard cosmogenic flux 
<Fq in the absence of BH production. For detection by the Auger ground array, Ces ~ 3.0, 


assuming maximal neutrino mixing and the (3 r value given above [|29|. The fluorescence 
detectors of HiRes provide additional sensitivity |3D|, as do those of Auger |53|. We conser¬ 
vatively take Ces = 3 for the combined rate in 5 years expected in the SM. Note, however, 
that the rate may be greatly suppressed for large BH cross sections, as anticipated. 

In contrast to Eq. ([H]), the rate for quasi-horizontal showers follows simply from Eq. 
and has the form 


A^qh — C< 


QH 


; cc 


+ cr 


BH 


4*0 CTqq 


(38) 


where Cqh = 2.5 for the Auger ground array, as noted previously. 

Given a flux and BH cross section crg H , both Aes and Aqh are determined. Event 
contours are given in the left panel of Fig. [12|. As can be seen, given a quasi-horizontal 
event rate Aqh, it is impossible to differentiate between an enhancement from large BH 
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FIG. 12: Left: Contours of constant number of quasi-horizontal showers Aqh (dashed) and Earth- 
skimming neutrino events Aes (dotted) as functions of source flux and BH production cross 
section Ug H . 5 year running times for Auger and HiRes are assumed. Right: Confidence level 
contours, assuming <& u = 4>q and <jg H = 3 <7q C , corresponding to (IVqh, Aes) ~ (10,0.2). 


cross section and large flux. However, in the region where significant event rates are ex¬ 
pected, the Aqh and Aes contours are more or less orthogonal, and provide complementary 
information. With measurements of Aqh and Aes, both cr^ n and may be determined 
independently, and neutrino interactions beyond the SM may be unambiguously identified. 
(See also Ref. 0.) 

As an example, consider the case in which a v GC = 3, and /<&$ = 1. On average, one 
would then observe a total of Aqh = 10 deep quasi-horizontal showers, an excess of 8 above 
SM expectations. On average, one also expects Aes ~ 0.2 Earth-skimming events. A SM 
explanation (with crg H = 0) of the deeply penetrating event rate would require = 4 

and predict 12 Earth-skimming events, a possibility that would be clearly excluded at high 
confidence level. 

More generally, one might try to salvage a SM explanation by attributing the observed 
rates to statistical fluctuations in both Aqh and Aes- Using a maximum likelihood method 
for Poisson-distributed data |86|, we give contours of constant \ 2 hi the right panel of Fig. [12. 
The possibility of a SM interpretation along the a^ H = 0 axis would be excluded at greater 
than 99.9% CL for any assumed flux. The power of the Earth-skimming information is 
such that the best fit is in fact found for < <1 >q! We find, then, that if even an excess 
of a handful of quasi-horizontal events is observed, by comparing to the Earth-skimming 
neutrino rate, attempts to explain the excess by SM interactions alone may be excluded. 
These arguments require only counting experiments, and do not rely on measurements of 
shower properties. 

BH production will most likely be accompanied by more model-independent sub- 
Planckian effects. In particular, neutral current neutrino cross sections may be enhanced 
in extra-dimensional scenarios through the exchange of KK gravitons. This will raise the 
quasi-horizontal rate, but will have very little effect on the Earth-skimming event rate, since 
neutrinos suffer very little energy loss during this process |L|. We expect such effects, then, 
to further enhance the ratio Aqh/Aes, making a SM explanation even more untenable. 
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So far, we have not explicitly considered the question of distinguishing BH events from 
other types of new physics. However, the prediction of enhanced quasi-horizontal event rates 
and diminished Earth-skimming rates is incisive. For example, new physics that increases 
quasi-horizontal rates by enhancing cross sections for uN —> tX will also increase Earth- 
skimming rates. The prediction of suppressed Earth-skimming rates relies on the efficient 
conversion of neutrino energy directly to hadronic energy, that is, a process with large 
cross section and large inelasticity. This is a peculiar property of BHs that separates BH 
production from other possible forms of new physics. The comparison between deep quasi¬ 
horizontal shower and Earth-skimming neutrino rates therefore not only effectively excludes 
a SM interpretation of BH events, but goes a long way toward excluding other new physics 
explanations. 


IX. SUMMARY OF RESULTS AND CONCLUSIONS 


In this work we have shown that cosmic ray observations in the recent past (AGASA) and 
in the near future (Auger) provide extremely sensitive probes of low-scale gravity and extra 
dimensions. We have focused on the production of TeV-scale BHs resulting from collisions 
of ultra-high energy cosmic neutrinos in the Earth’s atmosphere, and have considered the 
impact of various theoretical issues in the determination of the BH production cross section. 
In particular, mass shedding, the production of BHs with non-zero angular momentum, and 
a possible enhancement of the BH cross section can be expected to give minor perturba¬ 
tions. The exponential suppression proposed by Voloshin is more significant, but large and 
observable BH event rates are still possible. 

More specifically, in the case of n extra spatial dimensions compactified on an n-torus 
with a common radius, we have found the following: 


Present bounds on atmospheric BH production imply 95% CL lower limits on the 
fundamental Planck mass of Mb > 1.3 — 1.5 TeV for n = 4, rising to Mb > 1.6 — 
1.8 TeV for n = 7. These bounds follow from the non-observation of a significant 
excess of deep, quasi-horizontal showers in 1710.5 days of running recently reported 
by the AGASA Collaboration p5| . 


The absence of a deeply-penetrating signal in the Fly’s Eye data |2U[ also implies 
lower bounds on Mb- These are consistently weaker, however. For example, for 
n = 6, x m i n = 1, and the same (PJ) flux we have used, Ringwald and Tu find Mb > 
900 GeV [[HJ. We find this difference to be significant: the AGASA and Fly’s Eye 
constraints rely on identical theoretical assumptions, and given the scaling in Eq. (0. 
a factor of 2 difference in Mb bounds corresponds to a factor of more than 4 in 
acceptance or, equivalently, running time. 


The AGASA limits derived here exceed the D0 bound Mb ~ 0.6 — 1.2 TeV, where 
the variation reflects uncertainty from the choice of ultraviolet cutoff for graviton 
momenta transverse to the brane. The cosmic ray limits are subject to a separate set 
of uncertainties, discussed at length above, but follow from conservative evaluations of 
the neutrino flux and experimental aperture, and x min = 1. For x min = 3, these limits 
are somewhat reduced, but still generally exceed the Tevatron bounds. 


The cosmic ray bounds from AGASA therefore represent the best existing limits on 
the scale of TeV-gravity for n > 4 extra spatial dimensions. A summary of the most 
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FIG. 13: Bounds on the fundamental Planck scale Md from tests of Newton’s law on sub-millimeter 
scales, bounds on supernova cooling and neutron star heating, dielectron and diphoton production 
at the Tevatron, and non-observation of BH production at AGASA. Future limits from the Auger 
ground array, assuming 5 years of data and no excess above the SM neutrino background, are 
also shown. The range in Tevatron bounds corresponds to the range of brane softening parameter 
A /Mb = 0.5 — 1. The range in cosmic ray bounds is for x m \ n = 1 — 3. See text for discussion. 


stringent present bounds on Mjj for n> 2 extra dimensions is given in Fig. [H| 

• The reach of AGASA will be extended significantly by the Auger Observatory. If 
no quasi-horizontal extended air shower events are observed in 5 years (beyond the 
expected two SM neutrino events supplemented by as many as 10 hadronic background 
events), Auger will set a limit of Mp > 3 TeV, at 95% CL, for n > 4. Even in the case 
where the cross section is decreased by the exponential suppression factor in Eq. ©, 
a bound Mp > 2 TeV may be found under the same background assumptions. 

• Conversely, given the large reach of Auger, tens of BH events may be observed per year. 
We have discussed in some detail how combined measurements of quasi-horizontal air 
showers and Earth-skimming v T —> r events may be used to identify new neutrino 
interactions beyond the SM, even with complete uncertainty about the incident neu¬ 
trino flux. In the case of BH production, the quasi-horizontal event rate is enhanced, 
while the Earth-skimming rate is suppressed, since BH production in the Earth acts 
as an absorptive channel, depleting the SM rate. With counting experiments alone, 
one can therefore exclude a SM interpretation of BH events, and may distinguish BH 
events from almost all other possible forms of new physics. 

In conclusion, in the next several years prior to the analysis of data from the LHC, 
super-Planckian BH production from cosmic rays provides a promising probe of extra di¬ 
mensions. Searches for BH-initiated quasi-horizontal showers in the Earth’s atmosphere at 
AGASA provide the most stringent bounds on low-scale gravity at present, and the Auger 
Observatory will extend this sensitivity to fundamental Planck scales well above the TeV 
scale. 
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